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^ 1 Introduction 

These lectures are centered around a specific problem, the effect of weak repul- 

sive interactions on the transition temperature Tc of a Bose gas. This problem 
^ provides indeed a beautiful illustration of many of the techniques which have 

I been discussed at this school on effective theories and renormalization group. 

Effective theories ar(^ us(^d first in order to obtain a simple hamiltonian describ- 
ing the atomic interactions: because the typical atomic interaction potentials 
are short range, and the systems that we consider are dilute, these potentials 
can be replaced by a contact interaction whose strength is determined by the 
^ s-wave scattering length. Effective theories are used next in order to obtain 

I a simple formula for the shift in T^. one exploits there the fact that near Tc 

^ the physics is dominated by low momentum modes whose dynamics is most 

Q economically described in terms of classical fields; the ingredients needed to 

Q calculate the shift of Tc can be obtained from this classical field theory. Finally 

the renormalization group is used both to obtain a qualitative understanding, 
I and also as a non perturbative tool to evaluate quantitatively the shift in Tc- 

^ In the first lecture, I recall some known aspects of Bose-Einstein con- 

densation of the ideal gas. Then I turn to an elementary discussion of the 
interaction effects, and introduce an effective theory with a contact interac- 
tion tuned to reproduce the scattering length of the atom-atom interaction. I 
show that at the mean field level, weak repulsive interactions produce no shift 
in Tc- Finally, I briefiy explain why approaching the transition from the low 
temperature phase is delicate and may lead to erroneous conclusions. 
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In the second lecture I establish the general formula for the shift of Tc. 

H 

Cu where a is the s-wave scattering length for the atom-atom interaction, and 

Tc the transition temperature of the ideal gas at density n. This formula 
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holds in leading order in the parameter am}^^ which measures the diluteness 
of the system. Getting formula ([l]) involves a number of steps. First I explain 
why perturbation theory cannot be used to calculate AT^, however small a 
is. Then I show that the problem with the perturbative expansion is localized 
in a particular subset of Feynman diagrams that are conveniently resummed 
by an effective theory of a classical 3-dimensional field. The outcome of the 
analysis is the formula ([T]) where c is given by the following integral 

where the proportionality coefficient, not written here, is a known numerical 
factor, and S is the self-energy of the classical field, whose calculation requires 
non perturbative techniques. 

In the last lecture, I use the non perturbative renormalization group 
(NPRG) in order to estimate c. This requires the knowledge of the 2-point 
function of the effective 3-dimensional field theory for all momenta, and in 
particular in the cross-over between the critical region of low momenta and 
perturbative region of high momenta. This cross-over region is where the dom- 
inant contribution to the integral ([2| comes from. In order to obtain an accu- 
rate determination of S{p), it has been necessary to develop new techniques 
to solve the NPRG equations. Describing those techniques in detail would 
take us too far. I shall only present in this lecture the material that can help 
the student not familiar with the NPRG to understand how it works, and 
how it can be used. I shall do so by discussing several simple cases that at 
the same time provide indications on the approximation scheme that has been 
developed in order to calculate S{p). I shall end by reporting and discussing 
the results obtained with the NPRG, and compare them to those obtained 
using other non perturbative techniques. 

Recent discussion of Bose-Einstein condensation can be found in [Tl [51 [31 
m [S] . The equation ([ijfor the shift of Tc is derived and discussed in the series 
of papers [3 [H l^llU]. Much of the material of the last lecture is borrowed 
from the papers [HI US [13 [II [IS] . 



2 LECTURE 1. Bose-Einstein condensation 

2.1 Bose-Einstein condensation for the non interacting gas 

• The discussion of Bose-Einstein condensation of the ideal Bose gas in the 
grand canonical ensemble is standard. We consider a homogeneous system 
of identical spinless bosons of mass m, at temperature T . The occupation 
factor of a single particle state with momentum p is 



1 



(3) 
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where /i is the chemical potential. The number density of non-condensed 
particles is given by 

-^np = n{ii,T). (4) 

For small density, the chemical potential is negative and large in absolute 
value, e-'^/^ > 1. The gas is then described by Boltzmann statistics: 



.p ^ e -(4-M)/T^ ^ ^ e''/^A"3, (5) 
where A is the thermal wavelength: 



Unless specified otherwise, we use units such that fi = l,kB = 1- Boltz- 
mann's statistics applies as long as nX^ <C 1, that is, as long as the ther- 
mal wavelength is small compared to the interparticle distance. As one 
lowers the temperature, keeping the density fixed, the chemical potential 
increases, and so does the thermal wavelength. Eventually, as /i ^ 0_, the 
density of non condensed particles reaches a maximum 

d'p 1 . „ C(3/2) 

where Ci^) is the Riemann zeta-function and A the thermal wavelength. 
As one keeps lowering the temperature, particles start to accumulate in 
the lowest energy single particle state. This is the onset of Bose-Eisntein 
condensation, which takes place then when 

nX^ = C(3/2) w 2.612. (8) 

At this point the thermal wavelength has become comparable to the in- 
terparticle spacing. Eq. ([s]) defines the critical line in the n, T plane (see 
Fig.jlJ. In particular, the critical temperature is given by 

It is a function only of the mass of the atom, and of the density. For 
T < T", the chemical potential stays equal to zero and the single particle 
state p = is macroscopically occupied, with density uq oc n given by 

' rp \ 3/2N 



no(r) = n ( 1 - ( ^ ) ) . (10) 
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Fig. 1. Critical line in the {fJ.,T) plane for the ideal Bose gas 

The occurrence of condensation relies on the existence of a maximum of 
the integral Q fixing the number of non condensed atoms. This depends 
on the number of spatial dimensions. In 2 dimensions, the integral diverges 
logarithmically as /i — s- 0_: 



In this case, there is no limit to the number of thermal particles, hence no 
condensation in the lowest energy mode. 

The Bose-Einstein condensation of the non-interacting gas has some patho- 
logical features. In particular the compressibility 



diverges when /i ^ 0_. Also the fluctutation of the number of particles in 
the condensate, 



is of the same order of magnitude as the average value Nq. Such features 
can be related to the large degeneracy of states in Fock space that appears 
as /i — + 0. 




(11) 




(12) 



yZ(7V^ = ^NoiNo + 1) ~ A^t 



(13) 



Non Perturbative Renormalization Group and Bose-Einstein Condensation 



5 



Consider indeed the ground state in Fock space of the hamiltonian of the 
non interacting Bose gas, 



p 



(14) 



If < 0, then the ground state is the vacuum, with no particle present 
in the system: adding a particle costs a positive "free-energy" —fi. When 
H = 0, there appears a huge degeneracy in Fock space: all the states with 
an arbitrary number of particles in the single particle state p = are 
degenerate. If /x > 0, then there is no minimum, the more particle one 
adds in the state p = 0, the more free energy one gains {—fi for each 
added particle). Note that here the chemical potential does not fully plays 
its role of controlling the density: either it is negative, and the density is 
zero, or it is positive and the density is infinite; if it vanishes the density 
is arbitrary. 



To relate the large degeneracy to the large fluctuations (13 1, one may use 



a simple maximum entropy technique to determine the most likely state 
with A'o particles in the degenerate space. This calculation parallels the 
corresponding calculation done in the grand canonical ensemble at finite 
temperature. Maximizing 



5" = - ^ p„ In p„ + ^ ^ npn - 



Pn 



(15) 



with the constraints 



^Pn = 1, X! " 



(16) 



one finds 



For large A'o 



Z = 



/i = - In 1 



1 



-n/No 
Nn ' 



(17) 



(18) 



which leads indeed to the large fluctuations ( 13 1 



Of course, some of the pathologies of the ideal Bose gas in the grand canon- 
ical ensemble could be eliminated by working in the canonical ensemble, 
where the particle number is fixed (see e.g. |16p. However, we shall keep 
the discussion within the grand canonical ensemble, as it is closer to fa- 
miliar field theoretical techniques. This allows us in particular to treat 
Bose-Eisntcin condensation (of the interacting gas) as a symmetry break- 
ing phenomenon. 
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The large degeneracy of states in Fock space implies that infinitesimal 
interactions could have a large effect, and indeed they do: as soon as weak 
repulsive interactions are present the large fluctuations are damped, and 
the compressibility becomes finite . 



2.2 Interactions in the dilute gas 

• As we shall discuss later in this lecture, the dominant effect of the in- 
teractions in the dilute gas can be accounted for by an effective contact 
interaction whose strength g is proportional to the s-wave scattering length 
a: 

Zl 77-/7 

y(ri-r2)-5'5(ri-r2), g= . (19) 

m 

• In the mean field (Hartree-Fock) approximation, which is also the leading 
order in a, the effect of the interaction is a simple shift of the single particle 
energies: 

el^el + 2gn, (20) 
where the factor 2 comes from the exchange term. 

It is easy to see that this produces no shift in Tc- because the shift of the 
single particle energy is constant, independent of the momentum, eq. Q 
which gives the number density of non condensed particles, can be written 

n = n{i^- Aii,T) (21) 

with 

= 2gn. (22) 

Bose-Einstein condensation now occurs when fi — Aii = 0, instead of fi — 
in the non interacting case, but clearly the critical line is identical to that 
given by eq. (Isl. 



Of course, because the Hartree-Fock self-energy in eq. (20 1 depends on the 
density, the relation between the chemical potential at the transition and 
the critical density is more complicated than in the non-interacting case. 
The transition now takes place at a finite value of the chemical poten- 
tial, and some of the pathologies of the ideal gas disappear. In particular 
the compressibility remains finite at (fluctuations remain important 
though): 



dn _ T IpiT-pil + np) ^ 1 
dn 1 + ^ J np{l + np) ~ 2g 



2, r . . - ^- (23) 



• The presence of the interactions allows us to treat the phenomenon of Bose 
Einstein condensation as a symmetry breaking phenomenon. Let us return 
first to the ground state, and calculate the free energy 
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"^^-.n^ + lnl (24) 
whe re V is the volume. f2{n()) is the expe ctat ion value of the hamiltonian 



(14 1 to which is added the interaction (19 1, V ^ {g/2)aQaQaoaQ, in a 
coherent state containing an average number of particles iVo = noV in the 
state p = 0, |no) ~ exp ^/Noal\0) . In contrast to what happened for the 
ideal gas, it is now possible to minimize i7(no) w.r.t. uq in order to find 
the optimum ground state for a given /i. One gets 

no = fi/g. (25) 

Thus the ground state has now a finite number of particles, and the fiuc- 
tuations are normal, proportional to the square root of the mean value. 
The density in the ground state is n = dQ/d^ — hq, and the pressure is 
P = — f2{no = n/ g)/V = .gn^/2, so that also at zero temperature the com- 
pressibility is finite, = ndP/dn — gn? . Note however that the coherent 
state is a state with a non definite number of particles: the symmetry re- 
lated to particle number conservation is spontaneously broken. We shall 
return later to the similar picture at finite temperature, and come back to 
this issue of symmetry breaking. 



2.3 Atoms in a trap 

Although our main discussion concerns homogeneous systems, it is instruc- 
tive to contrast the situation in homogeneous systems to what happens in a 
trap. We shall consider here a spherical harmonic trap, corresponding to the 
following external potential 

V{r) = ^mw^r^. (26) 

• Consider first the non interacting gas. We assume the validity of a semi- 
classical approximation allowing us to express the particle density as the 
following phase space integral: 

"^"■^ " / (2^3 e(e(r,p)~t^)/T _ I ' (27) 

where e{r,p) = p^/2m+(l/2)ma;§r^. This requires that the temperature is 
large compared to the level spacing, fc^T 3> hcooi ^ condition well satisfied 
near the transition if iV, the number of particles in the trap, is large enough. 
The number density in the trap can be written as 

N = J d^rn{^i- ^mLu^r^,T), (28) 
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where T) is the function Q. Thus in a wide trap for which the semi- 
classical approximation is valid, the particles experience the same con- 
ditions as in a uniform system with a local effective chemical potential 



/i — imtjQr^. In particular, eq. (27 1 shows that the density at the center 
of the trap is related to the chemical potential by the same relation 
as in the homogeneous gas. It follows that T° is determined by the same 
condition as for the homogeneous gas, that is, n(0)A^ = 2.612 where n(0) 
is the density at the center of the trap. An explicit calculation gives 

(29) 



with C(3) — 1.202. Note that while the condensation condition is "univer- 
sal" when expressed in terms of the density at the center of the trap, the 
dependence of on iV depends on the form of the confining potential. 
In the present case, the N dependence can be obtained from the follow- 
ing heuristic argument. For a temperature T ^ the particle density is 
approximately given by the classical formula 



V2T 



(30) 



so that the thermal particles occupy a cloud of radius Rth ^ \JT jraui^ ~ 
a\g/X, where aho = ^/ y/mjjjQ is the characteristic radius of the har- 
monic trap and A ^ l/V mT is the thermal wavelength. The average 
density in the thermal cloud is n ^ NR^j^. At the transition, the in- 
terparticle distance is of the order of the thermal wavelength, so that, 
f7,i/3 ^ A^^/^A/a^Q ~ A"-^, or A^/^A^ ~ a\^, from which the relation 
r° ~ ujN^/^ follows. 

Note that the confining potential makes condensation easier than in the 
uniform case. This is related to the fact that the density of single particle 
states in a trap decreases more rapidly with decreasing energy than in a 
uniform system: it goes as e"*"^ in a trap and as e'^/^-i in a uniform system, 
where d is the number of spatial dimensions. It follows in particular that, in 
a trap, condensation can occur in d = 2, in contrast to the homogeneous 
case; the heuristic argument presented above yields then T° ~ uiN^^^. 
Note however that the effects of the interactions in a 2-dimensional trap 
are subtle (for a recent discussion see [T7] ) . 

The leading effect of repulsive interactions in a trap is to push the par- 
ticles away from the center of the trap, thereby decreasing the central 
density. This effect is analogous to that produced by an increase of the 
temperature. One expects therefore the interactions to lead to a decrease 
of the transition temperature. To estimate this, we note that, in the pres- 



ence of interactions, the density is still given by eq. (30 1 after substituting 
{l/2)mu!Qr'^ — > (l/2)ma;§r^ -|- 2gn{r). A simple calculation gives then the 
density n{r) (for r not too large, and to leading order in a) in terms of the 
density n{0) at the center of the trap [TO] : 
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n(r) « „(o)e-^("-o'-V2)(i-25„(o)/T) 
1 mujnr 



n{0) 



1 



2T 1 + 4aA2n(0) 



(31) 



This result suggests that the effect of the interaction can also be viewed 
as a modification of the oscillator frequency, Wq ujQ/{l + 4aX'^n{0). This 
is enough to estimate the shift in Tc: 

# = - - -aXMO) - (32) 

where, in the last step, we have used the fact that 77(0) ^ A^'^ at the 
transition. A more explicit calculation yields the proportionality coefficient 
-1.32 [18 . 

It is important to keep in mind that this effect of mean field interactions 
on Tc is very different from the one that leads to eq. ([T|) (indeed the sign 
of the effect is different). If we were comparing systems at fixed central 
density rather than at fixed particle number, there would be no shift of Tc 
(see the discussion in [W\). 



2.4 The two-body problem 

We now come back to the construction of the effective interaction that can be 
used in many-body calculations of the dilute gas. We shall in particular recall 
how effective field theory can be used to relate the effective interaction to the 
low energy scattering data. More on the use of effective theories can be found 
in the lecture by T. Schafer in this volume [19]. A pedagogical introduction is 
given in Ref. PU] . 

Consider two atoms of mass m interacting with the central two-body po- 
tential V{r), with r the relative coordinate and r — |r|. The relative wave 
function satisfies the Schrodinger equation 



i;{r) = Ei;{r). (33) 



The scattering wave function is given by {E = Ti^k'^ /m) 

7^k(r)=e''^-- + /(k',k) — , (34) 



r 



where k is the initial relative momentum (with k = |k|), while k' is the final 
relative momentum. The scattering is elastic and because the potential is ro- 
tationally invariant, the scattering amplitude /(k',k) is a function only of the 
scattering angle 9 in the center of mass frame, and of the energy E = h^k"^ /m. 
For short-range interactions, the interaction takes place predominantly in the 
s-wave, and the scattering amplitude becomes of a function of the energy only. 
It has then the following low momentum expansion: 
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/(^) = _l + mfc2_^fc = ^^E- (35) 

where a is the scattering length, r\) the effective range, and the neglected terms 
in the denominator involve higher powers of k. In the very low momentum 
limit, when fcro ka <C 1, one can ignore the effective range. Then the scattering 
amplitude depends on a single parameter, a. 

The scattering amplitude can be expressed as the matrtix element between 
plane wave states of the T-matrix: 

(k'|T|k) = -— /(k',k), (36) 

and the T-matrix itself can be calculated in terms of the Green function Go 

T = V- VGqT, (37) 

where 

Go{E) = (38) 

and E ^ E + ir] with E real for the retarded Green's function Gq. The 
following formal relations are useful 

H = Ho + V, G-^ = H- E = Gq^ + V, (39) 

and 

T = V -VGV -VGaT T-'^ = V-^+Ga. (40) 

Note that both G{E) and T{E) are analytic functions of E, with poles on the 
negative real axis corresponding to the energies of the bound states, and a cut 
on the real positive axis. 

Let us now turn to the many-body problem. Assuming that the atoms in- 
teract only through the two body potential V{r) , we can write the interaction 
hamiltonian as 

V=ll dh^dh^ V/(ri)V'^(r2)y(ri - r2)V^(r2)^(ri), (41) 

We demand that the local effective theory 

^-fS = % j «^'rV^(r)^t(r)^(r)V>(r), (42) 

reproduce the same scattering data in the two particle channel at low momen- 
tum as the original potential V{vx — v-i). We know that in the long wavelength 
limit the scattering amplitude depends only on the scattering length a, so we 
expect go to be related to a. 
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To establish this relation we calculate the scattering amplitude for the 
effective theory. For a contact interaction T(E) is given by 

^ ^ G^iE), (43) 



TiE) go 
where 

^°^^^ = /(S- pV^-^-..^ ^''^ 

To calculate Gq{E), which is divergent, we introduce a cut-off A on the mo- 
mentum integral 



— Em 27r^ 7o ^ ^ '^-^ ^ 



mA m^E f+°" , 1 



dp-^ E^— • (45) 



2' 9R=~Zr^ (46) 



It is then convenient to define a "renormalized" strength gji by 

1 1 TOyl 47ra 

— = 1" ' 9r= — 

9r .90 27r^ m 

where in the last relation a is the scattering length. This relation between gji 
and a is obtained by comparing T{E) obtained from eq. (44 1 {E > Q) 

1 1 m I — — 

= V i — VniE, (47) 



TiE) gn 4tt 



with eq. (35 1, and using the relation (36 1 to relate T{E) and f{E). 

Remark. One may improve the description by including the eflective range 
correction. This is done by adding to the hamiltonian a term of the form [23] 



9o 
2 



d^rV{iP'^ii})-V{iP'^iP), (48) 



and adjusting so as to reproduce the scattering amplitude of the original two 
body problem, at the precision of the effective range. At tree level in the effective 
theory, the calculation of the scattering amplitude yields 

(k3k4|rlki k2) = -2g'a [(ki - ka) ■ (ks - k4) + (ki - k4) ■ (ks - ka)] 

= 8(7^,fc^ (49) 

where in the second line k is the magnitude of the relative momentum. Note that 
in order to make the identification with the two-body problem discussed above, 
we have to pay attention that the two-body problem traditionally treats the two 
particles as distinguishable, whereas the present calculation involves matrix ele- 
ment of the T-matrix between symmetric two p arti cle states. Thus the tree level 
calculation of the T-matrix for the hamiltonian (41 1 yields {k3k4|r|ki k2) = 2go 
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which diflers by a factor 2 with the conventional definition of the scattering 
length. Staying with the usual convention, we therefore write 

y ^ 47r a ^ a ( ^ ^ aro , 2 

ml — roak^/2 + ika 

from which the identification of q'q follows 



ffo=5o-^- (51) 

For recent discussions on the application of effective field theory techniques 
to the Bose gas, see [HI [21 [22] • 



2.5 One-loop calculation 

Having at our disposal an effective many-body hamiltonian, we may now 
perform detailed calculations of the effect of the interactions. The one loop 
calculation that we present here gives us the opportunity to come back to 
the issue of symmetry breaking, illustrates the use of the delta potential and 
points to the difficulty of approaching the phase transition from below. 

The grand canonical partition function can be written as a path integral: 

Z = Tie-^"-''^ = [ V{4>,iP*) e-5 , (52) 



with 



(53) 

and P = l/T is the inverse temperature. The complex field ip{x) = 7/;(r, r) to 
be integrated over is a periodic function of the imaginary time r, with period 
p. The action S is invariant under a U{1) symmetry: 

V'(x) ^ e*"V(2^), Hx)* e-"'il;*{x). (54) 

It is convenient to add an external source j{x) linearly coupled to the 
bosonic field (thereby breaking the U{1) symmetry): 

Z[j]= f I?(^^V')e-5^ (55) 

where 

S^S = S- [ d^x [j{x)'ip*{x) + j*{x)ij{x)] . (56) 
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The loop expansion is an expansion around the field configurations which 
make the action stationnary. These field configurations, which we denote by 
4>q{x) are determined by the equation 



d_ 
dr 



2m 



<i>o{x) + 9a\<i>o?(t>o{x) = j{x) 



(57) 



and its complex conjugate. For uniform, time-independent, external sources 
j, eq. (57 1 admits constant field solutions (compare eq. (24 1): 

[-ti + go\M^]^o^j- (58) 

There are two solutions. For j — 0, (po — or |(/)op — fJ./go- The solution 
(^0 = corresponds to the vacuum state and is the stable solution when 
< 0. For fi > (and gQ > 0) the second solution corresponds to broken 
U{1) symmetry and BE condensation; in that case the solution 0o = is a 
maximum of the free energy. 

We now expand the field around the solution (f>o which corresponds to a 
minimum of the (classical) free energy: 



V'(x) — > (t>Q + iljix) il;{x)* 



4)a + tlj{x)* , 



(59) 



with t/'(a;) having only k ^ Fourier components, and keep in iS terms which 
are at most quadratic in the fluctuations. That is, we write Sj = 5o + 5i, with 
iSo the classical action 



So ^ 13V 



Ml0o|'+?(|0o|')'+J0S+r'/'o 



(60) 



and Si the one-loop correction 

r0 



Si 



dr / d^x 



ij*{x) 



d_ 



2m 



90 



{{^'^f ^x)^x) + {<iyo)' ^Hx)^Hx)) 



(61) 



The gaussian integral is standard and, after a Legendre transform to eliminate 
J, it yields the following expression of the thermodynamic potential: 



n 

V 



gnnl 



2V 



/3V 



where Ek is the Bogoliubov quasiparticle energy: 

El = (4 + 2gno - tif - {gn^f . 



(62) 



(63) 



Note that in this calculation we have used the relation ( 46 1 in order to replace 
the bare coupling constant by the renormalized one. This is the origin of the 
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fourth term in the r.h.s. of eq. (62 1 which ehminates the divergence in the 
sum over the zero point energies. Note that this replacement assumes that g^i 
and go are perturbatively related, which can only be possible if the cut-off A 
in (46l is not too large, i.e., if mAgu/ <C 1, or ^4 <C 7r/(2a). 

This calculation was first made in the context of Bose-Einstein con- 
densation by Toyoda [24], and led him to the erroneous conclusion that 
It is essentially a mean field calculation, the mean field 
being entirely due to the condensed particles. This approximation, equivalent 
to the lowest order Bogoliubov theory (see e.g. 125|), describes correctly the 
ground state at zero temperature and its elementary excitations. However, its 
extension near the critical temperature meets several difficulties; in particular 
it predicts a first order phase transition [26J, a point apparently overlooked 
in Ref. [24 . In fact, above T^, 4>o = 0, and the present one-loop calculation 
yields the free energy of an ideal gas, with no shift in the critical temperature. 



3 LECTURE 2. The formula for AT^ 

Let us start by considering again the phase diagram in Fig. 1. We want to 
determine the change of the critical line due to weak repulsive interactions. 
For small and positive values of the s-wave scattering length a, we expect the 
change illsutrated in Fig.[2j with the two critical lines close to each other. One 
can then relate, in leading order, the shift of the critical temperature at fixed 
density, AT^, to that of the critical density at fixed temperature, AtIc- Since, 
for a = 0, n° (X {T°f/'^, we have 

AT, 2 An, 

It turns out that it is easier to calculate at fixed temperature than at fixed 
density, and in the following we shall set up the calculation of Aric. 

In the previous lecture, we have seen that, under appropriate conditions, 
the effective hamiltonian is of the form 

H = Jd'r |v/(r) (J^ V(r) + |^t(r)^t(r)^(r)V>(r) 

(65) 

where go is related to the scattering length by (see eq. (46)): 



1 1 mA 47ra 

— = — + ^^ = • (66) 

9R 90 27r^ m 

This effective hamiltonian provides an accurate description of phenomena 
where the dominant degrees of freedom have long wavelength, kro <C 1, and 
the system is dilute, an^/'^ <C 1. Recall also that we shall be working in the 
vicinity of the transition where n^^^X ~ 1, with A the thermal wavelength. 
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Fig. 2. Expected change in the critical line caused by weak repulsive interactions 
characterized by a small scattering length a 



3.1 Condensation condition and critical density 

In order to exploit standard techniques of quantum field theory and many- 
body physics (see e.g.[25l|27l[2Hl[29]), we shall first relate the particle density 
to the single particle propagator. The particle density can be written as 

n = (■(/'^(r)V'(r)), 
= lim (TV'(T,r)V'^(0,r)> = lim G(t,0), (67) 

where 

V'^(T,r) =e^^V'^(r)e-^-^, ^{T,r) ^ e^" iP{r) e-''" , (68) 

are respectively the creation and annihilation operators in the Heisenberg 
representation and T denotes (imaginary) time ordering (see e.g. [H]). The 
single particle propagator is 

G(ri-T2,ri-r2) = (T^(ri, ri)7/.t(r2, ra)) 

= Trf— TV'(ri,ri)V'^(T2,r2)j. (69) 

It is a periodic function of r = ri — T2: for < r < /3, G(r — /3,r) — G(t, r), 
where P — l/T is the inverse temperature. Because of its periodicity it can 



16 Jean-Paul Blaizot 

be represented by a Fourier series 

G(t,p) = T^e-^""^G«(iu;„,p), (70) 

n 

where the a;„'s are called the Matsubara frequencies: 

Un = 2n7TT, (71) 

and we have also taken a Fourier transform with respect to the spatial co- 
ordinates. We are making here an abuse of notation: we denote by the same 
symbol the function and its Fourier transform, with the implicit understand- 
ing that the arguments, whether space-time or energy-momentum variables, 
are enough to specify which function one is considering. The inverse transform 
is given by 

G{iu;n,p)= / dre^""^G(r,p). (72) 
Jo 

In the absence of interactions, the hamiltonian is of the form 

Ho-iiN = J d^r ^t(r) (^-^ - V'(r) 

= ^ (e° - n) = ^ £p a^ap, (73) 



where 

p p 

In these formulae, V is the volume of the system, and the creation and anni- 
hilation operators satisfy [ap,ap,] = (5p,p'. The free single particle propagator 
can be obtained by a direct calculation. It reads 

Go^{iu)n, p) = £p - iu)n, (75) 

or, in imaginary time, 

Go(ti - T2,p) = (Tap(ri)a^(T2))o 

= e-^-'^^i-^^) [(1 + np)0(ri - T2) + np0(r2 - n)] , (76) 

where: 

rpj.g-/3(ffo-MA')cjta 1 

rip = {alap)o = ^^^_ff(H„_^jv'i " = ^Pe, _ 1 ' £p = £p - (77) 
Thus, for non interacting particles. 
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lim Go(t,p) = Up, (78) 

r — >U- 

so that the formula (67 1 yields the familiar formula Q of the density. 

The full propagator is related to the bare propagator by Dyson's equation 

G'"^(iw„,p) = G(7\iw„,p) + S{iujn,p), (79) 



where S(iuj,p) is the self-energy. In this case eq. (67 1 yields the following 
expression for the density 

n= lim ryf^ r^^^. (80) 

or equivalently the occupation factor in the interacting system 

np= lim T^e-*'^"^G(iw„,p). (81) 

n 

We shall approach the condensation from the high temperature phase. 
Then the system remains in the normal state all the way down to Tc- The 
Bose-Einstein condensation occurs when the chemical potential reaches a value 
such that (see e.g. [50]'): 

G-\lu = 0,p = 0)=0 or S{uj = 0,p = 0) = ^^. (82) 

At that point, 

G'\iiu„,p)=iLUn-el-[SiiiUn,p)-SiO,0)]. (83) 



By using the general relation ( 80 1 between the Green function and the 
density, one can then write the following formula for the shift Auc in the 
critical density caused by the interaction: 



Auc = lim T Ve-*""^ 



d^p 

(27r)3 



el + Siiuj„,p) ~ S{0, 0) - iu;„ - icj. 



(84) 



The first term in this expression is the critical density of the interacting system 
at temperature T, the second term is the critical density of the non interacting 
system at the same temperature. This formula makes it obvious that Auc 
vanishes if the self-energy is independent of energy and momentum. This is 
the case in particular when the interactions are treated at the mean field level, 
as we already observed. In this case, S — 2gn, and the formula above yields 

^"^ = / (^{e/5(eo+J-M_i-^]5?rrT} = '^' (^5) 

which vanishes since fic — '2.gn. 
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Fig. 3. The contribution to tiie self-energy E that is of leading order in a 



3.2 Breakdown of perturbation theory 

Because the interactions are weak, one may imagine calculating Arte by per- 
turbation theory. However the perturbative expansion for a critical theory 
does not exist for any fixed dimension c? < 4; infrared divergences prevent a 
complete calculation, as we shall recall. If one introduces an infrared cutoff fee 
to regulate the momentum integrals, one finds that perturbation theory breaks 
down when kc ^ a/A^, all terms being then of the same order of magnitude. 

The leading order on a is given by the diagram in Fig. [3] As we just saw, 
the contribution of this diagram to S is just the mean field value 2gn, and the 
net effect on Aric is zero. We shall then examine the second order contribution, 
given by the diagram displayed in Fig. |4] We shall see that this diagram is 
infrared divergent. Next, we shall show, using simple power counting, that 
such infrared divergences occur in higher orders and signal a breakdown of 
perturbation theory as one approaches the critical point. 



Second order perturbation theory 

The second order self-energy diagram is the lowest order diagram that is mo- 
mentum dependent and can therefore yield corrections to the critical density. 
It is displayed in Fig. |4] Its contribution is given by 



n'n" ^ ' ^ ' 

1 1 



(86) 

Anticipating that infrared divergence can occur, we focus on the contribution 
where ujn' — ^n" — 0, and calculate the difference S(iujn = 0,p) — S{iujn — 
0,p = 0), which we write from now on simply as Sci{v) — ^ci{Q)'- 

2^2 [ 



(27r)3 (27r)3 
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Fig. 4. The contribution to the self-energy S that is, a priori, of second order in a 



p+q 



(87) 



In this calculation we have used Hartree-Fock propagators, and set 
fi' = /i - Shf 



2m' 



where Shf — '^gn. The quantity fee, whose significance will appear more 
clearly as we progress in the discussion, plays the role of an infrared cutoff 
in the integrals. Note that fcc ^ (^' 0) when T T^. With this new 
notation 

el - // - {e + kl)/2m. (89) 

We shall also set 

U{p) = 2m{S,i{p) - S,i{Q)). (90) 
The integral in eq. ( [ST] ) can be calculated analytically (see [S]) and yields 

1 



arctan 

P ike 



■In 



1 



P 



- 1 



(91) 



This equation shows that U{p) is a monotonically increasing function of p, 
~ p2 at small p, and growing logarithmically at large p. This logarithmic 
behavior, obtained in perturbation theory, remains in general the dominant 
behavior of U{p) at large p, i.e., for fc^ ^ p^l/X. This equation (911 also 
reveals the anticipated infrared divergence: U{p) diverges logarithmically as 

The condensation condition (82), /i' = X!ci{0), reads 
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In 



A 



7 



(92) 



where the right hand side is —2m2Jci{0), 7 = 0.577... is Eulcr's constant, 
and the approximate equahty is vahd when A ^ kc, which we assume to be 
the case. Here A is an ultraviolet cut-ofF that we have introduced to handle 
the divergence of Sci (0) , and whose origin is the neglect of the non vanishing 
Matsubara frequencies. Typically, A ~ 1/A. Eq. (92) shows that kc ^ a/X^ 
when a is small 8 . 

This simple calculation illustrates the limits of a pertubative approach. 
The infrared cutoff kc aj}? introduces spurious a dependence. The conden- 
sation condition which relates the infrared cutoff to the microscopic length 
A, induces a spurious logarithmic correction which does not vanish as a — > 0. 
In fact such logarithms do appear as higher order corrections, as we shall see 
later, but axe absent in the leading order result. 



Higher orders 

The infrared divergences that we have identified in the second order calcu- 
lation persist, and worsen, in higher order contributions. This may be seen 
by using a simple power counting argument. Let us first consider diagrams 
in which all the internal lines carry vanishing Matsubara frequencies. We use 
again HF propagators so that all the functions that are integrated in the 
diagrams are products of fractions of the form 

{K^^e^-\ (93) 

where K denotes a generic combination of momenta; it is then natural to use 
the dimensionless products Kjkf. as new integration variables. Consider then 
a diagram of order a". The lowest order n = 2 has just been calculated, and, 
for large p/fc^, it is proportional to (a/A)^ ln(p//cc), where p is the external 
momentum. For n > 2, every additional order brings in one factor a from 
the vertex, one integration over three-momenta, a factor T, and two internal 
propagators. The contribution of the diagram can thus be written as: 

where -F is a dimensionless function, which we do not explicitly need here. 
The main point is that when one approaches the critical temperature, the 



coherence length becomes large so that the summation of terms (94) diverges. 
In the critical region, k^ ~ a/A^, so that all the terms in the perturbative 
expansion are of the same order of magnitude. Therefore, at the critical point, 
perturbation theory is not valid. 

Let us now assume that in a given diagram some propagators carry non- 
zero Matsubara frequencies so that one momentum integration (fc) will be 
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altered. For that integration, the presence of a non vanishing Matsubara fre- 
quency in the denominators of the propagators ensures that no singularity 
at fc = can take place. Essentially, in the corresponding propagators, fee is 



replaced by a term proportional to 1/A, so that one factor a/A fee in (94 1 is 
now replaced by a/ A. Compared to the diagram with only vanishing Matsub- 
ara frequencies, this diagram is down by a factor a/ A, and thus negligible in a 
leading order calculation of S. This reasoning generalizes trivially to diagrams 
containing more non vanishing Matsubara frequencies. 

Formula for Aric 

It follows from the previous discussion that in order to obtain the leading 



order shift in the critical density, one may retain in eq. ( 80 ) the contribution 
of the zero Matsubara frequency only. That is. 

An / f 1 1\ 

71" A^ Jo p^ + U (p) 



with U{p) given by eq. (90 1. Note that this integral is finite: U{p) ^ Inp at 



large p, and p + U{p) ^ p'^ at small p. Note also that since U{p) > (in 



fact the general qualitative behavior of U{p) is correctly given by eq. (91 1) 
the correction Arte is negative, implying a positive shift of Tc- 



3.3 Classical field approximation 

Once restricted to their zero Matsubara frequency components, the fields ip 
and V'^ can be considered as classical fields, and the entire calculation can be 
cast in terms of a classical field theory. To see that, let us expand the field 



variables in the path integral (52 1 in terms of their Fourier components: 



V;(r,r)=V'o(r) + T^e--"-V„(r), 



(96) 



where the w„'s are the Matsubara frequencies, and we have called "00(1") the 
r-independent part of the field. 



The partition function ( 52 1 can then be written as 



where ipo depends only on spatial coordinates, and 



(97) 
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Fig. 5. Some diagrams contributing to the parameters of the effective theory. The 
external legs represent the field ipo carrying vanishing Matsubara frenqency. The 
internal lines carry non vanishing Matsubara frequencies tJ„ . The first diagram from 
the left is a correction to the chemical potential; it is of order a, but does not con- 
tribute to the shift in T^. The second diagram is a correction to the i/'o coupling 
constant. The last diagram is the contribution to the ip^ vertex. The last two con- 
tributions are of order and respectively, and can be ignored in a leading order 
analysis. 



exp{~S,s[^Po]} = ^^J 'D^n^oe^^ (98) 



where S is the action (53). A/i and M2 are (infinite) normaHzation constants. 
The quantity S'offiV'o] is the effective action for tpQ. Aside from the direct 
classical field contribution to which we shall return shortly, this effective action 
receives also contributions which, diagrammatically, correspond to connected 
diagrams whose external lines are associated to ipo, and the internal lines are 
the propagators of the non-static modes ipn- A few examples are displayed in 
Fig. [sj Thus, a priori, i^efi ["00] contains operators of arbitrarily high order in 
ipQ. In the present case, however, it is easy to verify that the contributions 
beyond those kept in the classical action, are of higher order in a, and can 
therefore be ignored in a leading order calculation. 

The strategy which consists of integrating out the non-static modes in 
perturbation theory in order to obtain an effective three-dimensional theory 
for the soft static modes, is referred to, in another context, as "dimensional 
reduction" (see e.g. [5ni5^ [55]V 



In leading order, the effective action is the restriction of eq. ( 53 1 to the 



static mode ipo and the partition function (97 1 can be written as: 



A/" J V^Po exp j-Z^y" (H(r) - fiAf{r)) \ , (99) 



where V'o(r) is a three-dimensional field, and 

We shall refer to this limit as the classical field approximation. The zero Mat- 
subara component of the density is given by (|V'o(^)P)i which diverges in the 
effective theory. However, recall that by assumption, the wavenumbers of the 
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classical field are limited to k less than an ultraviolet cutoff yl ^ A ^ . In fact 
we shall not need to use a cut-off since the variation of the critical density is 



a finite quantity (see eq. (110 1 below) 



Remark. 

Ignoring the time dependence of the fields is equivalent to retaining only the zero 
Matsubara frequency in their Fourier decomposition. Then the Fourier transform 
of the free propagator is simply: 

Go(p) = — . (101) 



This may be obtaine d d irectly from ( 70 \ and ( 75 1 keeping only the term with 
uj„ = 0, or from eq. (1991). The classical field approximation corresponds also to 



(102) 



the following approximation of the statistical factor (see eq. ( |81[ l) 

1 T 



pPep 



Both approximations make sense only for ep <^ T, implying rip ^ 1. In this limit, 
the energy per mode is oc SpUp ~ T, as expected from the classical equipartition 
theorem. This long wavelength limit can also be viewed as a high temperature 
limit (the time dependence of the field becomes indeed unimportant as /3 — > 0). 
One should not confuse this classical field approximation with the classical limit 
reached when the thermal wavelength of the particles becomes small compared 
to their average separation distances. In this limit, the occupation of the single 
particle states becomes small, and the statistical factors can be approximated 
by their Boltzmann form: 

^ , ^ e'^'^"-") « 1, (103) 

where we have used the fact that e"*^*^ is large in the classical limit. 

At this point it is easy to understand the origin of the breakdown of 
perturbations theory. The critical region is characterized by the fact that all 



the terms in the integrand of ( 100 1 become of the same order of magnitude. 



This occurs for k ^ kc, with k^ such that: 

(104) 

where {T/ij/)kf is the contribution to the density of the modes with k ~ fc^. 
From eq. (104) we see that kc a/X^. For k ~ kc perturbation theory in 



a makes no sense, and in fact all terms in the perturbative expansion are 
infrared divergent. For fee ^ fc <C A^^, perturbation theory is applicable. 

Remark. In a trap the efi^ect of critical fluctuations is subleading. To see that 
let us estimate the size of the critical region in a trap [31]. Recall that at the 
transition the thermal wavelength is A ~ ahoN~^^^ , where a^o ~ 1 / y^mujo is 
the cha racteristic size of the harmonic oscillator trap (see the discussion after 



eq. |30|), and the size Rd of the thermal cloud where most of the particles 



sit is Rci ~ auoN^^'^ ~ ATV^''''. According to eq. (104), the critical region is 
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reached when the chemical potential deviates from its critical value by an amount 
Si fee /2m ~ a'^/(mA^). In a trap the effective local potential is of the form 
fi — imtJo''^- Taking fi at its critical value, one finds that the particles will be 
in the critical region as long as r ^ Rcr with mLOgRcr ~ a^/(mA*), or Rcr ~ 
a{aho/\y ~ aiV^''^. Thus the relative size of the critical region is Rcr/Rci ~ 
o/A, and the ratio of the particles in the critical region to the total number 
of particles is ~ (a/A)^. Under such conditions, it can be shown that one can 
use perturbation theory to estimate the corrections due to the interactions to 
the relation Q in a trap; the resulting contributions to the shift of Tc are then 
subleading in a, as compared to the mean field effect discussed in the first lecture 

m- 

In view of the forthcoming discussions, it is convenient to rescale the field 
tpo and to parametrize it in terms of two real fields (fi,tp2'. "00 = V mT{ipi + 
if 2)- The partition function then reads 

Z = J Vipe-^, (105) 

where S{(p) — {H — fiN)/T is given by: : 



5(^)= \-[dM^)]' + i,r^'{^) + ^W'{^)T\d'x, (106) 



where tp — tpf + and 



a 



r = -2mTn, u = 967r^-^. (107) 

In eq. ( |106[ ) we have kept the dimension d of the spatial integration arbitrary 
for the convenience of forthcoming discussions. The single particle Green's 
function G{p) is related to the inverse two-point function r^'^\p) of the clas- 
sical field theory by 

2mrG-i(p) = r(2)(p), p2^c/(p) ==r(2)(p) -r(2)(o). (io8) 

As it stands this field theory suffers from UV divergences. These are absent 
in the original theory, the higher frequency modes providing a large momen- 
tum cutoff ^ 1/A. This cutoff may be restored when needed, but, as we show 
later, since the shift of the critical temperature is dominated by long distance 
properties it is independent of the precise cutoff procedure. 

We shall also find useful to consider the 0{N) symmetric generalization of 



the Euclidean action (106 1. The field ip{x) then has N real components, and, 
e.g., 



N 



and the shift in the critical density is given by 
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N 



1=1 



with 6ij / r'^'^^k) the connected two-point correlation function. 

The advantage of this generahzation is that it provides us with a tool, the 
large N expansion, which allows us to calculate at the critical point (for a 
recent review see e.g. [35]). 

Linear dependence of the density correction 

It is now easy to see the origin of the linear relation between Aric and a. 



Note first that the action (1061 contains a single dimensionfuU parameter, 
u, r being adjusted for any given u to be at criticality. In fact the effective 
three dimensional theory is ultraviolet divergent, so there is a priori another 
parameter, the ultraviolet cut-off A ~ 1/A. It follows then from dimensional 



analysis that U{p) defined in eq. (90) can be written as 



U{p = xu)^u^(j{x,u/A). (Ill) 

Now, the diagrams involved in the calculation of U are ultraviolet convergent, 
so that U is in fact independent of the cut-off A, and the infinite cut-off limit 
can be taken. Note however that the validity of the classical field approxima- 
tion requires that all momenta involved in the various integrations are small 
in comparison with A ~ or, in other words, that the integrands are neg- 
ligibly small for momenta k ~ A^^. Only then can we ignore the effects of 
non vanishing Matsubara frequencies and use for instance the approximate 



form of the statistical factors (1021. In other words, the infinite cut-off limit 
is meaningful only if letting the cut-off becoming bigger than A^^ does not 
affect the results. This implies that uX ^ a/X is sufficiently small. 

In the region of validity of the classical field approximation, that is, for 
small enough u, a{x,u/A) becomes a universal function cr(a;), independent of 



u, and Aric in eq. (95 1 takes the form 



2u f , cr(x) , „, 

An,^ 7^ ^ V (112) 

ttA^ J a;^ H- a(x) 

showing that the change in the critical density is indeed linear in a. 



Renormalization group argument 



The linearity of the relation between the shift in Tc and the scattering length 
can also be understood from a simple renormalization group analysis. Let us 
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introduce a large momentum cutoff A ~ f/A, and a dimensionless coupling 
constant g 

g^A''-\o,(^^y'\ (ff3) 

At Tc the two-point function in momentum space satisfies the renormalization 
group equation ^37j 

(^ll + ^^^^ Wg ~ '^^^O ^^'^ ^ ■ ^^^^^ 

This equation, together with dimensional analysis, implies that the two-point 
function has the general form 

r(2)(p,^,g)^p2^(5)F(p/fc,), (ff5) 

where fcc = fcc(ff) is a function of g which, on d imen sional grounds, is propor- 
tional to A, so that AdAkc = fee- The ansatz (ff5) for F^^^p, A, g) provides 
then a solution of eq. (114) if Z{g)) and kc{g) obey the equations 

din Z{g) r,{g) 



dg /3(5)' 
91nfcc(g) 1 



(116) 
(117) 



dg P{g) ■ 

Since I3{g) = -(4- d)5+ (iVH-8)5V487r2 -f ©(g^)^ ^(^) is of order g for small 




Fig. 6. The /3-function with the two fixed points at g = and g — g* . 



g in d < 4; similarly ri{g) = {N + 2)g'^ / {72{8Ti'^f) -f 0{g'^). Therefore 
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(118) 



to leading order Z(g) =1. The function kc{g) is then obtained by integrating 
eq. (fTTTl, 



exp 



» , (J_ 1 
^ U(5')^(4-%' 



In d = 3, 



kcig) = A 



99 



9-9 



(119) 



(120) 



where g* is the infrared fixed point (see Fig. [6]). The scale kc{g) plays a spe- 
cific role in the analysis as the crossover separating a universal long-distance 
regime, where 



(121) 



governed by the non-trivial zero, g* , of the /3-function, from a universal short 
distance regime governed by the Gaussian fixed point, g — 0, where 



^(2) (p) cx p2 kc{g) p <^ A . 



(122) 



However such a regime exists only if kc{g) ^ A, i.e., if there is an intermediate 
scale between the infrared and the microscopic scales; otherwise only the in- 
frared behavior can be observed. In a generic situation g is of order unity, and 
thus kc{g) is of order A, and the universal large momentum region is absent. 
Instead fcc(ff) ^ A implies that g be small. Since g ~ a/A ^ 1, see eq. (1131, 
this condition is satisfied in the present situation. 

It is then easy to show, repeating essentially the same analysis as before, 
that with this condition, ATc oc kc{g). We set p — xkc{g), and find for the 
integral in eq. (110) the general form 



1 



1 



d'^p 

{2ttY Vr(2)(/c) p2 



= kc{g) 



1 



Fix) 



- 1 



(123) 



the g dependence is entirely contained in kc{g), and for small g, kc{g) ~ u. 

Recall that both the perturbative large momentum region and the non- 
perturbative infrared region contribute to the integrand in eq. or equiv- 

alently in eq. (1121, and that the functions F{x) or (t{x) cannot be calculated 



using perturbation theory. 

The 1 /N expansion allows an explicit calculation, and yields, in leading 
order for the coefficient c in eq. ([l]) the value c = 2.3 [7]. However, the best 
numerical estimates for c are those which have been obtained using the lattice 
technique by two groups, with the results: c = 1.32 ± 0.02 [SS] and c — 
1.29 ± 0.05 |39]. The availabilty of these results has turned the calculation 
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of c into a testing ground for other non perturbative methods: expansion in 
1/N [71135], optimized perturbation theory [321131], resummed perturbative 
calculations to high loop orders [151 . Note that while the latter methods yield 
critical exponents with several significant digits, they predict c with only a 10% 
accuracy. This illustrates the difficulty of getting an accurate determination 
of c using (semi) analytical techniques. 

Remark. Tlie linear dependence in a of the shift of Tc holds only if a is small 
enough, as we have already indicated. When a is not small enough various correc- 
tions need to be taken into account that alter the simple linear law. In particular, 
corrections come from the non vanishing Matsubara frequencies, and their im- 
pact on the effective theory for i/jq. Such corrections have been analyzed in detail 
in [33] (see also [§]). The net result is the following expression for AT^. 

^ = c{an^''') + [c^ ln(an^''') -I- c^'] {an^^'^ f + 0{a'n). (124) 

Aside from such corrections, the evaluation of ATc in higher order requires that 
one improves the accuracy of the effective hamiltonian and include for instance 
effective range corrections. 



4 LECTURE 3 : The Non Perturbative Renormalization 
Group and the calculation of c 

The analysis of the previous lectures has shown that the coefficient c in the 
formula ([T|) can be written as: 



2567r-^ 





(C(3/2))^/^ 



(125) 



where ^{fi) represents the change, due to interactions, of the fluctuations of 
the scalar field of the 0{N) symmetric model, and is given by: 



Nu 



27r2 x2 + cr(x) ' 



(126) 



where a{x) — u~'^U{p = xu) and U{p) is essentially the self-energy of the field 
at criticality (see eqs. ( |110| ) and ( |108| )). In order to get the numerical factor 
in eq. (|125|), we have combined eqs. (|l]), ([64|), {[sf and (|107|). Eq. \l2h\ has 



been written for arbitray N in order to be able to compare with all available 
results. Bose-Einstein condensation corresponds to = 2. 



The integrand of eq. ( 126 1, calculated in the approximate scheme described 
in |llj , is shown in Fig. |7| The minimum occurs at the typical scale kc which 
separates the scaling region from the high momentum region where perturba- 
tion theory applies. As we have already emphasized, the difficulty in getting 



a precise evaluation of the integral (1261 is that it requires an accurate de- 



termination of a{x) in a large region of momenta including in particular the 
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-4e-07 



0.01 



- 0.0001 



integrand 



le-06 



le-08 



- le-10 



-5 

log(p/u) 



Fig. 7. The function a{x = p/u) at criticality, and the integrand of eq. ( 1261 as a 
function of ln(p/u), for u ^ lO^"* for which (t(x) is independent of it. From Ref. [llj . 



crossover region between two different physical regimes; as we have seen in the 
previous lecture, this cannot be done using perturbation theory, eventhough 
the coupling constant, ~ a, is small. 

In this last lecture, I shall show how the non perturbative renormalization 
group (NPRG) can be used to calclate a{x). The NPRG 01 US US |13 
|48] (sometimes called "exact" or "functional" , depending on which aspect 
of the formalism one wishes to emphasize) stands out as a promising tool, 
suggesting new approximation schemes which are not easily formulated in 
other, more conventional, approaches in field theory or many body physics. It 
has been applied successfully to a variety of problems, in condensed matter, 
particle or nuclear physics (for reviews, see e.g. [HJ [SHI EI]). In most of these 
problems however, the focus is on long wavelength modes and the solution 
of the NPRG equations involves generally a derivative expansion which only 
allows for the determination of the n-point functions and their derivatives 
essentially only at vanishing external momenta. This is not enough in the 
present situation where, as we have seen, a full knowledge of the momentum 
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dependence of the 2-point function is required. We have therefore developed 
new methods to solve the NPRG equations. The following sections will briefly 
present some of the ingredients involved, without going into the technical 
details of their implementation. These can be found in Refs. pTllT^IT^IHlllSj . 
from which much of the material presented here is borrowed. For definiteness, 
we focus on the 0{N) symmetric scalar field theory with action (1061. See 
also the lecture by H. Gies in this volume }52j for further introduction to 
these techniques, and their application to other theories, in particular non- 
abelian gauge theories. Other applications of the renormalization group to 
Bose-Einstein condensation can be found in [551IM1I55] . 



4.1 The NPRG equations 

The NPRG allows the construction of a set of effective actions Fi^lcl)] which 
interpolate between the classical action S and the full effective action In 

[<j)] the magnitude of long wavelength fluctuations of the fleld is controlled 
by an infrared regulator depending on a continuous parameter k which has the 
dimension of a momentum. The full effective action is obtained for the value 
K = 0, the situation with no infrared cut-off. In the other limit, corresponding 
to a value of k of the order of the microscopic scale A at which fluctuations 
are suppressed, /^k=a[<?!'] reduces to the classical action]^ 

In practice the fluctuations are controlled by adding to the classical action 



(1061 the 0{N) symmetric regulator 



^^.M ^M) R.{q) ^.(-9), (127) 

where Rk. denotes a family of "cut-off functions" depending on k. As we just 
said, the role of AS^ is to suppress the fluctuations with momenta q ^ n, 
while leaving unaffected those with k. Thus, typically, R^iq) when 
q <^ K, and i?K(9) — > when q^ n. There is a large freedom in the choice of 
Rniq), abundantly discussed in the literature [5^1 1571 [SSJ [5^. The choice of 
the cut-off function matters when approximations are done, as is the case in 
all situations of practical interest. Most of the results that will be presented 
here have been obtained with the cut-off function proposed in |58j : 

i?„(q)cx(K2_q2)g(^2_g2)^ (^28) 

This regulator has the advantage of allowing some calculations to be done 
analytically. 

For each value of k, one deflnes the generating functional of connected 
Green's functions 



^ Note that depending on the choice of the regulator, not all fluctuations may 
be suppressed when k = A. However, for renormalisable theories, and if A is 
large enough, the effects of these remnant fluctuations can be absorbed into a 
redefinition of the parameters of the classical action. 
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W^[J] ^IuJd^ exp - AS^[^] + J 



(129) 



The Feynman diagrams contributing to are those of ordinary perturba- 
tion theory, except that the propagators contain the infrared regulator. We 
also define the effective action, through a modified Legendre transform that 
includes the explicit subtraction of AS^' 



(130) 



rM = -w^iM + / d^'x 0(x) j^x) - ASM, 



where is obtained by inverting the relation 



SJ{x) ' 



(131) 



for (f)^j{x) — (p. Note that, in this inversion, </) is considered as a given variable, 
so that J(j, becomes implicitly dependent on n. Taking this in to account, it is 
easy to derive the following exact flow equation for /^k [0]: 



1 



d^R^iq) G{k, g), 

where the trace tr runs over the 0(iV) indices, and 

Gr.i(g; k; 0) = r^fi^; q- 0) + 5,^RM- 



(132) 



(133) 



with r^"^^ the second functional derivative of with respect to (f) (see eq. ( 134 ) 



below). Eq. (1321 is the master equation of the NPRG. Its solution yields the 
effective action = /^k=o[0] starting with the initial condition r'K=/i[0] = 
5 [(/)]. Its right hand side has the structure of a one loop integral, with one 
insertion of df^Rt^{q^) (see Fig. [8|. This simple structure should not hide the 
fact that eq. ( 132[ ) is an exact equation (G in the r.h.s; is the exact propagator), 
and as such it is of limited use unless some approximation is made. Before we 
turn to such approximation, let us further analyze the content of eq. ( 132[ ) in 
terms of the n-point functions. 

As well known (see e.g. [37]), the effective action is the generating 
functional of the one-particle irreducible n-point functions. This property ex- 
tends trivially to /^k [(/>]: 



E 



^ j d^x^ ... j d^xr, ri") [</.; Xi, • • • , x„] <P{x^) . . . 0(x„). 



(134) 



By differentiating eq. (132 1 with respect to </>, letting the field be constant. 



and taking a Fourier transform, one gets the flow equations for all n-point 
functions in a constant background field 0. For instance, the flow of the 2- 
point function in a constant external field reads: 
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Fig. 8. Diagrammatic illustration of the r.h.s. of the flow equation of the effective 



action, eq. (132 1. The crossed circle represents an insertion of OkRk, and the thick 



line a full propagator in an arbitrary background field. 



d, 



q,-p-q;K; 



(3) 

xGkiiq + p; 0)-fh,„(-p,p + q, -q; k; (j))Grni{.q\ n; 



1 



Gtj{q; k; (j))rj^tlkip, -p, q, -q; k; (l))Gkiiq; k; < 



(135) 



The corresponding diagrams contributing to the flow are shown in Fig. [9] 





Fig. 9. The two diagrams contributing to the flow of the 2-point function in a 
constant background field. The lines represent dressed propagators, Gk- The cross 
represents an insertion of di^Rk- The vertices denoted by black dots are ri^' and 

(4) 

. When the background field vanishes, so does the diagram on the left. 



When the external field vanishes, this equation simplifies greatly since then 
r^^^ vanishes, and F^"^^ is diagonal : 



(136) 



which defines the self-energy S (this U differs from S^i introduced in sect. 3.2 



by a factor 2m; see e.g. eqs. ( |90| ) and (108l). We get then 



d.rl^l\n;p) = -\j ^dMq)GH^;q)r!^tlM;p,-p,q,-q), (137) 



where 
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G-\n,q)=q^ + RM + S{^^;(l)- (138) 
Similarly, the flow of the the 4-point function in vanishing fleld reads: 



X [g{k; q + pi +P2)-r'ib]j(K;pi,P2, q, ■)ri%{n\p?„Pi, -g, ■) 

+ G{k; q + pi +p3)r^'^i^{K;pi,p3, q, ■)rl2j{K;P2,P4:, -q, ■) 
+ G{K;q + pi+ J34)/id]j (k; Pi, Pi, q, ■)r'^uj P3 , P2, -g, •) } 

(139) 

where we have used the abbreviated notation r^^lj{K;pi,p2,q,—pi — P2 — q) 



d'^q 



r^flj{K;pi,p2,q,-)- The four contributions in the r.h.s. of eq. (1391 are 
represented in the diagrams shown in figs. |10| and [TT] 



Eqs. (137) and ( |139 1 for the 2- and 4-point functions constitute the be- 



ginning of an infinite hierarchy of exact equations for the n-point functions, 
reminiscent of the Schwinger-Dyson hierarchy, with the flow equation for the 
n-point function involving all the m-point functions up to to = 2. Clearly, 
solving this hierarchy requires approximations. A most natural approxima- 
tion would rely on a truncation, in order to close the inflnite hierarchy, as 
commonly done for instance in solving the Schwinger Dyson equations. For 



instance, one could ignore in eq. ( 139 1 the effect of the 6-point function on the 
flow of the 4-point function, arguing for instance that the 6-point function is 
perturbatively of order . However such truncations prove to be not accurate 
enough for the present problem. Besides, the NPRG offers the possibility of 
exploring new approximation schemes that involve the entire functional r[(f>] 
rather than the individual n-point functions. An example of such an approxi- 
mation is provided by the derivative expansion described in the next section. 



4.2 The local potential approximation 

The derivative expansion exploits the fact that the regulator in the flow equa- 



tions (e.g. eqs. ( 137 1 or ( 139 1) forces the loop momentum q to be smaller than 
ft, i.e., only momenta q ^ k contribute to the flow. Besides, in general, the 
regulator insures that all vertices are smooth functions of momenta. Then, in 
the calculation of the n-point functions at vanishing external momenta pi, it 
is possible to expand the n-point functions in the r.h.s. of the flow equations 
in terms of q^ j , or equivalently in terms of the derivatives of the field. 

In leading order, this procedure reduces to the so-called local potential 
approximation (LPA), which assumes that the effective action has the form: 
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Fig. 10. Diagrammatic illustration of the r.h.s. of the flow equation for the 4-point 
function, eq. ( |139[ ): contribution of the 4-point functions (represented by black disks) 
in the three possible channels, from left to right. The crossed circle represents an 
insertion of OkRk, and the thick line a full propagator. 




Fig. 11. Diagrammatic illustration of the r.h.s. of the flow equation for the 4-point 
function, eq. ( 139 1: contribution of the 6-point function F^-^^ (represented by a black 
disk). The crossed circle represents an insertion of OkRk, and the thick line a full 
propagator. 



(140) 



where p = 4)i4>i/2. The derivative term here is simply the one appearing in the 
classical action, and Vf^{p) is the effective potential. The exact flow equation 
for Vk is easily obtained by assuming that the field (p is constant in eq. ( 132 1. 
It reads: 



1 r d'^q 

2 J (2^ 



= I TWZudMq) {(N - 1)Gt{k; q) + Gl(«; q)} , (141) 



where Gt and Gl are, respectively, the transverse and longitudinal compo- 
nents of the propagator: 

G.,(«;9) = GT{K;q) - + GU^;q)^. (142) 



By using the LPA effective action, eq. (1401, one gets 
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G^\K-q)=q^ + V'{p) + Rk{q), 

GI\k- q)=q^ + V'ip) + 2pV'\p) + RM. 



(143) 



with V'{p) = dV/dp and V"{p) = d^V/dp^. 

The solution of the LPA is well documented in the literature (see e.g. 
[SIT, '35'). In this section, we shall just, for illustrative purposes, solve approx- 
imately these equations, by keeping only a few terms in the expansion of the 
effective potential in powers of p (thereby effectively implementing a trun- 
cation which ignores the effect of higher n-point functions on the flow) . The 
derivatives of Vi^{p) with respect to p give the n-point functions at zero ex- 
ternal momenta as a function of k. We shall introduce a special notation for 
these n-point functions in vanishing external field: 



dp 



d'V. 



p=0 



dp2 



d^V^ 



p=0 



dp^ 



(144) 



p=0 



The equations for these n-point functions are obtained by differentiating 
once and twice eq. ( 141 1 with respect to p, then setting p = 0, and using the 



definitions in eq. (1441. One gets, respectively 

{N + 2) 



(2) 



and 



where we have defined 



(145) 



(146) 



d'^q 
(27r)° 



Kd^R^{q^)G^{n-q) 



the explicit form in the second line being obtained for the regulator ( 128[ ) and 



K- 



2d-i T.'in dr{d/2). 



(148) 



We have used the fact that for vanishing external field, the propagator is 
diagonal and S{K.,q) = ml (see Eqs. ( [l36| ) and ([T38|)). Eqs. ( [l45| and ([l46| 
are solved starting from the initial condition ai n = A: 



9A 



u 
3' 



(149) 



In order to factor out the large variations which arise when k varies from 
the microscopic scale A to the physical scale k = 0, and also to exhibit the 
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fixed point structure, it is convenient to isolate the explicit scale factors and 
to define dimensionless quantities: 

If one assumes that <C 1, and ignore the contribution of /ik, then the 
equation for becomes: 

K^^{d-4)g^ + 2iN + 8)gl, (151) 
and can be solved explicitly: 

ff« = fvi^d^ (152) 

where g* is the value of g at the IR fixed point, g* ~ (4 — d)/ {2{N + 8)), and 
Kc the value of k for which = .g*/2. We have: 



KcA'*-^ g*-g^ g* 



^ ' 9 A 9 A 



(153) 



where the last approximate equality holds if g* 3> cja- In this regime, one 
recovers the qualitative feature already discussed in the previous lecture: there 
exists a well defined scale Kc ^ A, k'^~'^ = uKd/{3g*), that separates the 
scaling region, dominated by the IR fixed point, where g = g* , from the 
perturbative region, dominated by the UV fixed point g = (when k 3> 
Kc, one can expand g^, in powers of Kc/k; in leading order g^ = (w/3)(l — 

The local potential approximation, or a refined version of it, enters in an 
essential way in the approximation scheme that we have developed in order 
to calculate S{p)- Further insight can be gained by studying the ri-point 
functions in the large N limit, to which we now turn. 



4.3 Correlation functions in the large AT limit 

The LPA, as well as the higher orders of the derivative expansion, give accurate 
results for the correlation functions and their derivatives only at zero external 
momenta. In order to get insight into the effect of non vanishing external 
momenta we consider now the correlation functions in the large N limit (at 

fixed uN) [601 EH Ea in]- 

For vanishing field, the inverse propagator has the same form as in the 
LPA, eq. ( |138[ ), where the running mass is here given by a gap equation 

Nil f d'^n 

= ^ + ^ y {G{n; q) G{A- q)). (154) 
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Fig. 12. Diagrams representing the contribution to r^2M{i^)PijP2,P3,P'i) in the 
first channel of eq. ( 155 1. This diagram represents also a typical contribution to the 
function gK.{p) of eq. (1561. External momenta are counted as incoming momenta, 
so that pi + p2 flows into the bubble chain, and pa + P4 = — (pi + P2)- 



The 4-point function has the following structure: 

Sl2h49K{Pl +P2) + ^13'5245k(P1 +P3) + '5i4'5235k(P1 + Pa) , 

(155) 



where S'k(p) is given by 

n (r,\ - 

with 



i?d(Ac;p)= J ^GiK;q)GiK;p + q). (157) 

Finally the 6-point function r^23immi'^'^ Pii P21 P3t PAt Qt "l) (summation over 
repeated indices is understood) is of the form 

A^234mm ('^ ; Pi , P2 , PS , P4, Q, -q) 

= hniPl +P2)5l25zi + hf,{pi +P3)(5l3<^24 + ^k(P1 + P4:)Sud23, (158) 

with 

K{p)^Ng,{0)gl{p) J ^^G\n-q')G{K-q' + p). (159) 

Note that a single function, gK,{p), suffices to calculate all the n-point functions 
in the large N limit. These can be obtained in a straightforward fashion by 



calculating the corresponding Feynman diagrams with a regulator (see Figs. 12 



and 13 1. It is however easy to verify that the various n-point functions that 
we have just written are indeed solutions of the flow equations in the large N 
limit. 
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Fig. 13. Diagrams representing the contribution to -ri234mm('^iPiiP2,P3,P4, g, — g) 
in eq. (158 1. The small circles represent the function g^ip) with p the momentum 
flowing through the vertex. The convention for momenta is as in Fig. |12[ 



To this aim, one notes first that eq. (137) reduces to an equation for the 
running mass: 



1 



= -2^5k(0) J ^dMQ)G'{K;q), 



(160) 



and using eq. (1561, it is easy to check that this equation is equivalent to the 
gap equation, eq. (154). 

Next, we observe that in the large TV limit, a single channel effectively 
contributes in eq. ( 139 ) for the 4-point function; one can then use the following 
identity in this limit: 



^1% {i^',PiiP2, q,-q-pi -P2)r^tlj ps ,P4,-q, q-pa -Pi) 

= NgliPl +P2)Sl2S34, 



together with eq. (158) for and one obtains: 

^d.g.ip) = Ngl{p)jf\n-p) - ^h^{p)lf\n), 



(161) 



(162) 



('2) 

where the function ' {n) is that defined in eq. (147 1, with n = 2 . The 

'd 

d'^q 



function {i^]p) is obtained from the general definition 



{2n) 



jKdMq)G''-\i^;q)G{K;p + q). 



(163) 



Note that J^J'\k;p = 0) = ^"^(k). 

At this point we remark that the flow equation for g^ip) can also be 
obtained directly from the explicit expression (156), in the form: 

N 



d'^q 



^k5k(p) = j j;^G{K;q)G{K;q + p) 



(164) 
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It is then straightforward to verify, using eqs. (1601 and (1591 that eqs. (162) 



and (164 1 are indeed equivalent. The first term in eq. (1621 comes from the 
derivative of the cut-off function in the propagators in eq. (1641, while the 



second term, which involves the 6-point vertex, comes from the derivative of 
the running mass in the propagators. 

Note that eqs. (1601 for to^ and (1621 for gi^{p = 0) become identical 
respectively to eqs. (145) and (146) of the LPA in the large N limit, a well 
know property [5^ . 

Let us now analyze characteristic features of the function (p) . For sim- 
plicity we specialize for the rest of this subsection to c? = 3. Furthermore, for 
the purpose of the present, qualitative, discussion, one may assume ttIk — 0. 



This allows us to obtain easily gnip) from eq. (156) in the two limiting cases 
p = and K = 0. In the first case, we have 



5k (0) 



1 



3 1 



uN 1 
97r2 K 



This is identical to eq. (152), with here g* 
the other case, we have 



l/(2iV) and 



7k=o(p) 



1 



P 



3 1 



uN 1 
48 p 



3p + Pc 



(165) 
Nu/9'k'^ . In 

(166) 



with Pc = uN/A8. 

One sees on eqs. (165) and (166) that the dependence on p of g^^oip) is 
quite similar to the dependence on k of g^ {p = 0). In particular both quantities 
vanish linearly as k ^ or p ^ 0, respectively. The result of the complete 
(numerical) calculation, including the effect of the running mass ( i.e., solving 
the gap equation ( |154[ ) and calculating gnip) from eq. (156l), can in fact be 
quite well represented (to within a few percents) for arbitrary p and k by the 
following approximate formula 



9k{p) 



X 



31 + X 



X 



K 



P_ 
Pc 



(167) 



This simple expression shows that p, when it is non vanishing, plays the 
same role as k as an infrared regulator. In particular, at fixed p, the flow 
of 5k (p) stops when X becomes independent of k, i.e., when k ^ p{kc/pc), 
with KcIpc = 16/37r^ ~ 0.54. This important property of decoupling of the 
short wavelength modes is illustrated in Fig. [14] As shown by this figure, and 



also by the expression (167), the momentum dependence of the 4-point func- 



tion can be obtained from its cut-off dependence at zero momentum. In fact 



Fig. 14 suggests that, to a very good approximation, there exists a parameter 
a such that g(K;p) « g(K;0) when k > ap, and g{K;p) ~ g{K = ap;0) when 
K < ap. From the discussion above, one expects a ~ k^/pc — 16/37r^ « 0.54, 
which is indeed in agreement with the analysis in Fig. |14| 
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Fig. 14. The function i7k(p) (in units of A) obtained from a complete numerical 
solution of eqs. (160 1 and (162 1, as a function of k/u (in a logarithmic scale) for 
five values of p: from bottom to top, p/u = 0.001,0.01,0.1, 1 and 10. The envelope 
corresponds to p = 0. This figure illustrates the decoupling of modes: for each value 
of p, the flow stops when k ^ ap. The various horizontal asymptotes (dotted Uines) 
correspond to the single value a = 0.54. From Ref. (12) . 



The decoupling of modes can also be visualized in Fig. [15] The p- 
dependence of J^"'*(k;p) (eq. (163)) shown in this figure is relatively simple: 
whenp <C K, J^\k;p) ~ I^^\k); whenp k, J^\k;p) vanishes as 1/p^. On 
a logarithmic scale the transition between these two regimes occurs rapidly at 
momentum p ~ k, as illustrated on Fig. 15 A similar analysis can be made 



for the contribtution of the 6-point function in eq. (1621, and one can show 



that the ratio of the p-dependence of the second term in the right hand side 
(the one proportional to /i^) is, whenever it is significant, proportional to that 
of the first term, the proportionality coefficient being a function of k only. 



All this suggests that one can rewrite eq. (162) for g^ip) as follows 



d.g.{p) « Ngl{p)0{\ 



9 9 

a p 



)/f (^)(1-F0, 



(168) 



where a is a parameter of order unity. The 0-function ensures that the flow 
exists only when k > ap, and stops for smaller values of k. 

These approximations, together with another one that we shall discuss 
more thoroughly in the next section have been used in order to construct 
approximate flow equations in Refs. [T2j[T3]. 
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4.4 Beyond the derivative expansion 

The arguments on which the derivative expansion is based can be generalized 
in order to set up a much more powerful approximation scheme [14j that we 
now briefly present. This scheme allows us to obtain the momentum depen- 
dence of the n-point function with a single approximation. 

We observe that: i) the momentum q circulating in the loop integral of 
a flow equation is limited by k; ii) the smoothness of the n-point functions 
allows us to make an expansion in powers of j k} , independently of the value 
of the external momenta p. Now, a typical n-point function entering a flow 

equation is of the from Pi"'' (pi , p2 , •■•,Pn-i + Q,Pn — <!'■, 0), where q is the loop 
momentum. The proposed approximation scheme, in its leading order, consists 
in neglecting the g-dependence of such vertex functions: 

ri"'Hpi,P2, ■■■,Pn-i + q,Pn~ q; 0) ^ r'i"^(pi,P2, ...,p„-i,p„; (169) 

Note that this approximation is a priori well justified. Indeed, when all the 



external momenta vanish pi = 0, eq. (1691 is the basis of the LPA which, as 



stated above, is a good approximation. When the external momenta pi begin 



to grow, the approximation in eq. (169) becomes better and better, and it is 



trivial when all momenta are much larger than k. With this approximation. 



eq. (1371 becomes (for simplicity we set here iV = 1) 
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-jG„(g2; 0)^(4) (p, 0, 0; (/.)G,(<z2; 0)1 . (170) 



Now comes the second ingredient of the approximation scheme, which 
exploits the advantage of working with a non vanishing background field: 
vertices evaluated at zero external momenta can be obtained as derivatives of 
vertex functions with a smaller number of legs: 



(n+l) 



(pi,P2, •■•,Pn,0; 



d(j) 



(171) 



By exploiting eq. ( 171 1, one easily transforms ec 



(recall that and Ik are related by eq. (133l): 



( 170 I into a closed equation 



d.ri^\p^-4>) = J ^,dMq')Gliq'; 
x<,^^^^^i%^VG«((P + .)^;0)- 



ld^d^\p, ~p; 



.(172) 



Note the similarity of this equation with eq. ( 141 1 : both are closed equa- 



tions, the vertices appearing in the r.h.s. being expressed as derivatives of the 
function in the l.h.s.. 

The approximation scheme presented here is similar to that used in |12lll3j . 
There also the momentum dependence of the vertices was neglected in the 
leading order. However further approximations were needed in order to close 
the hierarchy. The progress realized here is to bypass these extra approxima- 
tions by working in a constant background field. 

The resulting equations can be solved with a numerical effort comparable 
to that involved in solving the equations of the derivative expansion. The 
preliminary results obtained so far are encouraging |15j . 



4.5 Calculation of A{ip^) 

We come now to the conclusion of these lectures, where results concerning the 
numerical value of c obtained with the NPRG will be presented. 

The approximation developed in [TTJ [T5J [TH] is based on an iteration 
scheme, where one starts by building approximate equations for the n-point 
functions, that are then solved exactly. To give a flavor of this method, consider 
more specifically the equation ( 139 1 for the 4-point function. An approximate 
flow equation is obtained with the following three approximations: 
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Our first approximation is that discussed in sect. 4.4: we ignore the q 
dependence of the vertices in the flow equation, i.e., we set g = in the 
vertices r'^^^ and F^^^ and factor them out of the integral in the r.h.s. of 



eq. (139 1. The second approximation concerns the propagators in the flow 



equation, for which we make the replacements: 

G{p + q}^GLPA'{q)0(^l'^^ (173) 

where a is an adjustable parameter. A motivation for this approximation may 
be obtained from the analysis of the n-point function presented in sect. 4.3 



(see eq. (1681 and Fig. 15 1. This approximation introduces a dependence of 
the results on the value of a. The third approximation concerns the function 
r(6) for which we use an ansatz inspired by the expressions of the various 
n-point functions in the large N limit (see again the discussion at the e nd of 



sect. 4.3): one assumes that the contribution of /^^^^ to the r.h.s. of eq. (139) 
is proportional to that of the other terms, the proportionality coefficient F,^ 
being only a function of k. The same proportionality also holds in the LPA 
regime, which allows us to use the LPA to determine F^. 

These three approximations result then in the following approximate equa- 
tion for 

X |©(k^ - a^(pi +P2f)) r/*^„(K;pi,p2,0, -pi -p2) 

xF^fl^{K,;p3,p4, -P3 - P4, 0) + permutationsj . (174) 

This equation can be solved analytically in terms of the solution of the LPA 
(actually a refined version of it). This is done by steps, starting form the mo- 
mentum domain a'^{pi +p2)^,ct^(pi + Pa)^! ck^(pi +P4)^ ^ i^^-, where it can 
be verified that the solution is that of the LPA itself. The solution of this 



equation is then inserted in eq. ( 137 1 for the 2-point function and leads to the 
leading order determination of S{p). The scheme is improved through an iter- 
ation procedure. At next-to-leading order, one improves P*-'*-' by establishing 
an approximate equation for F''^\ The solution of this equation i s the n used 



in eq. (1391 for F'^^\ and the resulting F'^^'^ is used in turn in eq. (1371 to get 
S at next-to-leading order. 

The results obtained with this method for the value of c are reported in Ta- 



ble 1 and Fig. 16 for various values of N for which results have been obtained 
with other techniques, either the lattice technique [SHIESIES], or variationally 
improved 7-loops perturtbative calculations |42j (other optimized perturba- 
tive calculations have also been recently performed, and are in agreement with 
those quoted here; see [SJIID]). The results reported in Table 1 have been ob- 
tained with an improved next-to- leading order calculation where, among other 
things, the parameter a of eq. ( |173 1 has been fixed by a principle of minimum 



sensitivity (see fTS] for details). A recent (approximate) calculation along the 
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c 


N = 1 


iV = 2 


iV = 3 


iV = 4 


N = 10 


TV = 40 


N = oo 


lattice [38J 




1.32 ±0.02 












lattice |39] 




1.29 ±0.05 












lattice [03] 


1.09 ±0.09 






1.60±0.10 








7-loops [H] 


1.07 ±0.10 


1.27 ±0.10 


1.43± 0.11 


1.54± 0.11 








large A'' [3 














c = 2.33 


NPRG 


1.11 


1.30 


1.45 


1.57 


1.91 


2.12 





Table 1. Summary of available results for the coefficient c. The last line contains 
the results obtained in [l3] . 



lines described in sect. 4.4 yields, for iV = 1, the value c — 1.2 [15]. Let me 
also mention the result c = 1.23 obtained, for iV = 2 in Ref. [64]; however, as 
discussed in [T3] it is difficult to gauge the quality of the approximation made 
in [Ml- 

For < 4, where we can compare with other results, the values of c ob- 
tained with the present improved NLO calculation are in excellent agreement 
with those obtained from lattice and "7-loops" calculations. What happens at 
large values of N deserves a special discussion. As seen in Fig. [16] the curve 
showing the improved leading order results extrapolates when A'^ — > oo to a 
value that is about 4% below the known exact result [7]. A direct calculation 
at very large values of A'^ is difficult in the present approach for numerical rea- 
sons: since the coefficient c represents in effect an order 1/7V correction (see 
[J), it is necessary to insure the cancellation of the large, order A^, contribu- 
tions to the self-energy, in order to extract the value of c. This is numerically 



demanding when A'^ ^ 100. Fig. 16 also reveals an intriguing feature: there 



seems to be no natural way to reconcile the present results, and for this matter 
the results from lattice calculations or 7-loop calculations, with the calcula- 
tion of the 1/A'' correction presented in Ref. [36]: the dependence in 1/N of 
our results, be they obtained from the direct NLO or the improved NLO, 
appear to be incompatible with the slope predicted in Ref. [35] for the 1/A^ 
expansion. 
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